When is Wegner's flow geodesic? 
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Abstract. Wegner' s method of flow equations offers a powerful tool for diagonalizing 
a given Hamiltonian and is widely used in various branches of quantum physics. Here, 
generalizing this method, a condition is found, under which the corresponding flow of a 
quantum state becomes geodesic in a submanifold of the projective Hilbert space. This 
implies the geometric optimality of the present method as an algorithm of generating 
stationary states. The result is illustrated by analyzing some physical systems. 
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In recent years, much attention has been focused on Wegner's method of flow 
equations [1,2]. This method offers a powerful tool for diagonalizing a given quantum- 
mechanical Hamiltonian, which has been devised in analogy with the theory of 
renormalization groups and has widely been applied to a variety of problems in 
condensed matter physics, nuclear and particle physics, and quantum information. 
Examples include the effective Hamiltonian of the Anderson impurity model [3], a 
Dirac particle in an external electromagnetic field [4], an effective spin- spin coupling 
arising from spin-phonon chains [5], the t — t Hubbard model for high-temperature 
superconductivity [6], the Tomonaga-Luttinger model for interacting spinless electrons 
in one dimension [7], localized superfluidity [8], the Lipkin-Meshkov-Glick model in 
nuclear physics [9], quantum phase transition in the interacting boson model [10], light- 
cone quantum chromodynamics [11], electron-electron and electron-phonon interactions 
in the Hubbard- Holstein model [12], and stationary photon-atom entanglement [13]. 

Wegner's method employs a continuous unitary transformation represented by the 
operator, U{1), where / e [0, °°) is referred to as the flow parameter. U{1) transforms 
the original Hamiltonian H = H(0) to H(l) = U(l) HU\l), which satisfies 
d H(l) I dl H(l)], where r](l) is the anti-Hermitian generator given by 

ri(l) = [dU(l) / dl]U ] (I). Wegner's choice for r]{l) to diagonalize (or, block- 
diagonalize) the Hamiltonian reads 

7]{l) = T] w {l)^[H i {l),H .S)l (1) 
where H d (l) and H Q _ d (l) stand for the diagonal and off-diagonal parts of H{1), 
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respectively. It can be shown [1,2] that H _ d (l) tends to vanish in the limit /— >°°. 
Mathematically, this procedure corresponds to the Jacobi algorithm for eigenvalue 
problems, which defines a steepest-descent flow in the space of matrices with the 
Frobenius norm [14]. 

Although Wegner's method has been applied to many problems as mentioned above, 
its inherent properties do not seem to be explored well. It is our opinion that this method 
is more than just a mathematical tool for diagonalizing a Hamiltonian. The purpose of 
this paper is to reveal an interesting geometric property hidden behind Wegner's flow 
equations. For it, we first generalize Wegner's method itself. The original method of 
Wegner is a special case within this generalized framework. Next, we consider the flow 
of a given quantum state through a submanifold of the projective Hilbert space 
composed of rays explained below. Then, we find a condition, under which the flow 
becomes geodesic in the submanifold. This remarkable property is illustrated by 
analyzing some physical examples. Our discussion may also have significance for 
laboratory experiments, since one often need realize a stationary state from a certain 
initially prepared state via quantum operations [13]. 

Let us start our discussion with considering a quantum system in a J-dimensional 
Hilbert space, where d is either finite or infinite. Its unitary-transformed Hamiltonian, 
H(l), is decomposed into two parts: H(l) = H d (l) + H _ d (l), where H d (l) and 

H o-d(f) = Y^-i H o-a(0 (l^A<d-l) are the diagonal and off-diagonal parts, 
respectively. Using the complete set of the normalized eigenstates of H d (l), 
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i \u „ ) \ , we write them in the following forms : 

V n /} n =l,2,...,d 



»d(0 = £e.(/)|» 11 X»-|. 



(2) 



re=l 



CI- (/) 



(3) 



where e n (l) is the nth eigenvalue of H d (l) and (/)'s are the complex expansion 
coefficients. The nonzero index i a describes the off-diagonality and is ordered without 
loss of generality as follows: < j\ < i 2 < ■ ■■ < i A . It is understood that C ( n " + ) (0 = 
(C!, (/) = 0), if n + i a >d{n-i a <\). 

Here, we generalize Wegner's method as follows. Instead of taking the whole of 
H _ d (I), we employ only a single term, say H ( a) d (l): 



n {a \i) = [H A {i\ h ( :1(d]. 



(4) 



Clearly, Wegner's choice is r] w (/) = £ fl _y a) (0- It can be found after a straightforward 
calculation that the corresponding generalized flow equation, 



dHjl) 
dl 



= [r 1 ia, (l), H(l)], 



(5) 



leads to (d/dl)Y d 

^^m, n=l, m^n 



H(l) 



= -*LL [^ +ia (l)-e n (l)] 2 \c^ a (l) 



implying that the off-diagonal elements of HQ) tend to decay as / increases. Therefore, 
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the present method generalizes Wegner's in the sense that it uses only a single off- 
diagonal term. 

Next, let us translate the flow of the Hamiltonian into the flow of a state. Given a 
normalized state \lff), the stationary Schrodinger equation reads H\lff(°°)) = E\\i/(°°)), 
where |y(0) = U ' (l)\y/) . As well known, two states different from each other only by 
total phases are equivalent in quantum mechanics, and therefore a physical state is 
represented by a ray. Accordingly, the quantum-state space is the projective Hilbert 
space. Eq. (5) determines the flows of the physical coefficients contained in the 
Hamiltonian, which depend on the parameters appearing in the unitary operator. Let us 
explicitly write as follows: U=U(a), = U\a)\y), where a= (a 1 , a 2 , a k ). 

The set of parameters, a, defines a local coordinate on the submanifold of the 
projective Hilbert space, which is referred to as the quantum evolution submanifold [15]. 
Then, the Fubini-Study metric [15-17] induced on this submanifold is, up to the second- 
order infinitesimals, given by 

ds 2 = l-\(y/(a)\y/(a+ da))\ 2 =g ij (a)da l da J , (6) 

where the metric tensor is expressed in terms of the anti-Hermitian operator 
G,(a) = [a,[/(a)] U\a) (3,. =d/da l ; i = 1,2,. ..,k) as follows: g u (a) = 
-(U2){y\G i (a)Gj(a) + Gj (a)G i (a)\yf) + {y/\ G i (a) | y/)(y/ 1 Gj (a) | iff) . Here and 
hereafter, Einstein's convention is understood for the repeated indices. The equation for 
a geodesic curve parametrized by the arc length, a (s), is given by 
d 2 a h /ds 2 +T"j (da 1 1 ds)(da j /ds) = 0, where T*. is Christoffel's symbol defined 
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by r hij = 8hl r l ij = (i/2) (a ,. g hJ + a j g ih - d hgij ). 

Let us parametrize the curve by the flow parameter, instead of the arc length, i.e., 
a(l) , and consider the functional, S = 2 dl L with L(a, d) = ds I dl = ^gijd' a 1 

( d' = da' (I) I dl), which is the arc length in the interval [l v l 2 ]. The variation of S with 
respect to a (I) is calculated to be 8S I da' = -(1 / 4) L" 3 X t , where 



X : = 



+2 



XYl \dl da'\ 



da 1 



Y )-2( ¥ \%\ ¥ ){ Y \£L\ Y ) 



1, 



dr] 

M 1 



W)\{{¥\rfW)-{wHw) 2 ) (7) 



with ri=ri(a(l)) = [dU(a(l))/ dl]U ] (a(l)) = G i (a(l))a'(l), provided that the 
following relation has been used: {d I dl)[dr]l dd') = dr] I da' + \r], dri/da']. 

We wish to prove that under a certain condition the generalized flow equation (5) 



defines a geodesic curve associated with the "initial state" \y/) = \u n ), which is an 
eigenstate of H d (l). (This is the state of relevance, because we are considering the flow 
to an exact stationary state.) That is, we are going to show that, under the condition 
found later, the above X { vanishes for r\ = ri {a) in Eq. (4). 

Clearly, (u n |?7 <fl) and its derivatives with respect to / and d' vanish. Furthermore, 

it can be shown by using Eqs. (2)-(4) that (u n | r\ (a)2 , drj (a) I dd' j|i< n ^ also vanishes. 
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Therefore, X t for \y/) = \u n ) and rj = ri (a) is reduced to 



X, = 2lu M a)2 \u n )(u n 



dr] (a) dri (a) 



dl ' 3d' 



,(a)2 



\dn I w ,a?7 ( 

W„ 7— W„ )(K„ h— rl«, 



n / \ n \ 



3d' 1 " 



(8) 



This equation can be rewritten by using Eqs. (2)-(4) again as follows: 



X,. = 4 



+ |d (0) IT <5 + . 

I n | / v n,m+i a 



Re 



m=l 



n+i a n+l a . n n o 

3d 7- ~di wxThx n,m+i ° 



D 



dl 



f g 

j; t—t n,m+i a 

u 1 m=l 



D 



3d'' 3d 



^■1 f ,5 

m=l 



(9) 



where =(e B+j (/)- E n (l))Cl+i (/). The /-derivatives appearing in this expression 

should be calculated from the generalized flow equation in Eq. (5). It is also noticed that 

a' and d' are not independent in Eqs. (8) and (9) any more, since Eq. (4) has already 
been used. 



Now, we present our finding. If the label a satisfying 



i a .*2i a ,3i a (!<«'< A) 



(10) 



can be taken for rj (a) (l), then sandwiching the flow equation in Eq. (5) by three sets, 



u\ and 



and \u n ), (u 



and 



u n _: ), we find 



r7 (a) (/)X^ 1 ^ ) (0l"„_ ! >0, (u n+ia W a \l),^__H^(l) 



u„ ) = 0, and 



«„_, j = 0, respectively. Then, the sandwiched flow 
equations are (i) dC ( n a \l) I dl = -(e n (l) - £„_,. V))* C ( f(l) , (ii) dC ( n ^{l)ldl = 

-(£ n+lfl (0-e„(0) 2 C^,(/), and (iii) (e n+l (l) + £_,(/) -2£„(/))cr ) (0 C«(l) = 0, 

respectively. It is noticed that, from (i) and (ii), the phases of C ( n "\l) and C ( n "^ (I) are 
found to be independent of /. (Case-A): If C[ a \l) = . (/) = 0, then X,. in Eq. (9) 
obviously vanishes. (Case-B): If C ( n a) (l) = and C^. (/) * (or, C ( n a) (l)*0 and 
C*^. (/) = 0), then X t vanishes, due to the fact that the phases of D l f(l) and D^. (/) 
are independent of /. (Case-C): If both C { n a) (l) and C ( n a + ) (I) are nonzero, then (iii) 
yields £ n+i (/) + £„_, (/) - 2e n (l) = 0. Combining this relation with (i) and (ii), we 



obtain a crucial relation that 



C^l) (I) is proportional to I C\ a> {l) I. And, such a relation 



makes X ( vanish again. Therefore, we conclude that the flow equation with the 

generator rj {a) (l) with the label a satisfying the condition in Eq. (10) gives rise to the 
geodesic flow of \u n j. This is the main result of the present work. 

In what follows, we illustrate our result by analyzing some physical examples. 

The first example we consider is the generalized harmonic oscillator. The 
Hamiltonian reads H = (Oa ] a + Xa 1 + Xa 2 + V. a' and a are the creation and 
annihilation operators obeying the algebra: [a, a ' ] = 1, [a \ a f ] = [a, a] = 0. 0) (> 0) and 
V are real constants, whereas X is a complex constant. The relevant unitary operator for 

this system is the squeezing operator: U(l) = exp^(l)a fl - %*(l)a 2 j / 2j, where the 



complex coefficient §(/) is parametrized as £(/) = r(/) e~ 2 ' W) (0 < r(/), < 0(/) < 27r). 
Accordingly, (a 1 , a 2 ) = (r, 0). The condition U(0)-I (with the identity operator I) 
leads to r(0) = 0. The "initial state" is taken to be the number state, \n) = (nl)~ U2 a 1 "|0), 
with the normalized ground state satisfying a|0) = 0. The corresponding metric is [15]: 

ds 2 = (1 / 2) (n 2 + n + 1) [dr 2 + (sinh 2 2r)d(j) 2 ], which shows that the manifold is the 

Lobachevsky space. The transformed Hamiltonian is written as follows: 
H(l) = co(l)a 1 a + Ml) a f +X (l)a 2 + v(l) . The coefficients appearing here depend not 
only on their original values at / = but also on £(/). In this case, 
H _ d (l) = Ml)a l +X(l)a 2 is the one and only off-diagonal part. Therefore, the 

generator in Eq. (4) is identical to Wegner's choice, T] w (l) = 2o)(l) A(/)a + 2 - A*(/)a 2 J, 

and the condition in Eq. (10) is automatically fulfilled. The flow equation for is 
given by d X (I) I d I = -4ft) 2 (I) X (I), showing that the phase of does not depend on 
/. Then, one can explicitly find that Eq. (8) vanishes. Comparing [dU(l) I dl] t/ f (/) 
with ri w (l), we obtain (j)(l) = const, which in fact turns out to make both 8SI 8r{l) 
and 8S/8(j)(l) vanish, where S is the arc-length functional defined in terms of the 
above metric. We also mention that the spectrum of H d (l) = (0(l)a ] a + v(l) is equally 
spaced, and accordingly Case-C is realized. Thus, Wegner's flow is geodesic. 

The second example is a spin-*, S = (S x , S y , S z ), in a constant external magnetic field, 

B. The Hamiltonian reads H = S ■ B in an appropriate unit. The unitary operator to be 
considered is U(l) = exp [a (/) S + - (J* (I) where S ± =S x ±iS y and the complex 
coefficient a{l) is parametrized as oil) = [0(1) 1 2]e~ mi) (0 < 9(1) <7i,0< (p(l) < 2k) 
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with 0(0) = 0. The basic commutation relations satisfied by the spin operators are as 
follows: [s„, S+] = ±S ± , S_] = 2S z . The local coordinate is given by 

(a 1 , a 2 ) = (6, q>). The initial state is taken to be \m) s (m = -s, - s + 1, 0, s - 1, s), 
which satisfies S z \m) s = m \m) s . The metric is found to be given by [15]: 

ds =(l/2)(s +s-m 2 )[dd 2 + (sin 2 d)d(p 2 \ which is of a sphere. The transformed 

Hamiltonian is written as H(l) = f3 z (l)S z + f3(l)S + + f3*(l)S_. f5 z (l) is real, whereas 
[3(1) is complex. They depend not only on B but also on (7(1). The one and only off- 
diagonal part is H _ d (l) = fi(l)S + + f3*(l)S_, and so Wegner's choice is employed for 
the generator: i] w (l) = p z (/)[/* (/) S + -/T(/)S_]. Therefore, clearly the condition in Eq. 
(10) is fulfilled. The flow equation for j8(Z) is d (5(1) I dl = -f5 2 (l)f5(l), from which the 
phase of fi(l) is seen to be independent of /, and accordingly Eq. (8) vanishes. 
Comparison of the above generator, r\ w (l), with [d U(l) I dl] U ' (/) yields (p(l) - const, 
which leads to the fact that the variations of the arc-length functional (S, calculated 
using the above metric) with respect to 6(1) and (p(l) vanish. Also, the spectrum of 
H d (l) = f3,(l)S z with s > 1 is equally spaced, and so realized is Case-C. Thus, Wegner's 
flow is geodesic (i.e., the great circle). 

The third and last example is the Jaynes-Cummings model [18], which describes a 
two-level atom interacting with a single-mode radiation field. The Hamiltonian is given 
by H = a> a 3 / 2 + a>a l a + K(a + a + a_a l ). Here, co , (O, and k are the transition 
frequency of the atom, the frequency of the radiation, and the coupling constant, 
respectively, a 1 and a are the creation and annihilation operators of the radiation field 
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satisfying the same algebra as in the first example, ct's are the operators of the atom, 
which are given in terms of the orthonormal basis, the ground state \g) and the excited 
state \e), as follows: o + = \e)(g\, o _ = \g)(e\, <7 3 = \e)(e\-\g)(g\. This model is known 
to be exactly solvable. The associated unitary operator is U(l) 

= LI=o[ a » (/)<T ^-l n X n l + Pn-i(f)<T-<r + \n)(n\ + y n (l)o + \n)(n + l\ + 8 n (l)a_\n + l)(n\], 

where \n) = (n!) _1/2 a tB |0) is the n-photon state (for the details of this unitary operator, 

l l 2 l l 2 

see Ref. [13]). The unitarity of this operator leads to the conditions: CM + y I =1, 

|a B | = |£»l' W = W and a nK + P* n Y„=0- Fram the P olar forms > «„ = 
|a n |exp^'0 a j and so on, follows the condition on the phases: a +0 p —B y -9 S 

= {2m + X)K (m - 0, ± 1, ± 2, ...). Here, we choose \e)\n) as the initial state. |a„|, a , 
6^ , and Q y form as a set of independent coordinate variables. The metric is ds 2 - 

I i 2 \ _1 I I 2 I I 2 / I \ 2 \( \ 2 
1- a I I d\cc n \ +\(X„\ \1-\cxj \[dd a -d0 Y I . The 0^ -dependence disappears 

due to the above choice of the initial state. The transformed Hamiltonian is H{1) - 
£~ =o [A n (I) o + o_ | n){n \ + B n _, (I) o_o + \ n){n \ + C n (I) a + \ n){n + 1 1 + C* n (I) a_\n + l)(n |] , and 

so H _ d (l) = ^" =Q [C n (l)c + \n)(n + l\ + C* n (l)c_\n + l)(n\], where the coefficients, A„(/)'s 

and 5 n _!(/)'s are real, and C n (/)'s turn out to be also real, as can be seen from the flow 
equations. They are expressed in terms of the physical coefficients contained in the 
original Hamiltonian as well as the coefficients appearing in U (I) (for details, see Ref. 

[13]). Let us consider Wegner's choice: J] w (l) = £~_ Q [A B (/)- B n (l)] C n (l) 
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x (<7 + |n)(n + l|-(7_|n + l)(n|]. Then, we find that Eq. (8) vanishes. The flow equations 
give rise to the result, dO a (1)1 dl = dG y (l)ldl = 0, which explicitly makes the 

variations of the arc-length functional with respect to | «„(/) , a (/), and 9 y (I) all 
vanish. In addition, comparing r\(l) = [H d (l), H _ d (l)] with r\ w (I) above, we see that 
Case-B is realized. Thus, again Wegner's flow is geodesic. 

In conclusion, generalizing Wegner's method of flow equations, we have found the 
condition, under which the corresponding flow of a quantum state becomes geodesic in 
a quantum evolution submanifold. We have illustrated this by employing the 
generalized harmonic oscillator, the spin in an external magnetic field, and the Jaynes- 
Cummings model. The present result implies that the method of flow equations is not 
just a mathematical tool for diagonalizing a Hamiltonian but provides an optimal 
strategy in quantum state engineering for realizing a stationary state from a given initial 
state. 
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